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O ■ Abstract 
m 

Bodeker has recently argued that non-perturbative processes in very high 
■ temperature non-Abelian plasmas (such as electroweak baryon number vi- 

(— i , olation in the very hot early Universe) can be quantitatively described, to 

Qh| leading logarithmic accuracy, by a simple diffusive effective theory. Bodeker's 

Q_ij effective theory is intended to describe the long-distance transverse electric 

and magnetic fields which are responsible for non-perturbative dynamics. His 
effective theory, however, also contains long-wavelength longitudinal electric 
fields. We discuss several subtleties in the treatment of longitudinal dynamics 
which were not closely examined in Bodeker's original treatment. Somewhat 
to our surprise, we find that within its domain of validity Bodeker's effective 
theory does correctly describe both longitudinal and transverse fluctuations. 
We also show that, as far as the transverse dynamics of interest is concerned, 
Bodeker's effective theory could be replaced by a transverse-only theory that 
removes the longitudinal dynamics altogether. In the process, we discuss 
several interesting aspects of stochastic field theories. 



I. INTRODUCTION 

Non-perturbative processes in a hot non-Abelian plasma at or near equilibrium are as- 
sociated with slow evolution of magnetic gauge fields.^ The characteristic spatial scale R of 



1 This is explicitly argued in ref. l]], but this fact is also implicit in earlier analysis of specific thermal effects 



such as plasmon damping rates o: 



fast-moving particles and the color conductivity ||. 



non-perturbative gauge field fluctuations and the associated time scale t for their evolution 
are of order 



1 1 . 

R ~ -r^;, t ~ — 7-——, — — , (1.1a) 

for small coupling. Alternatively, the characteristic spatial momentum k and frequency uj 
are 

k~g 2 T, u~ g *Tln(l/g). (1.1b) 

For a review, see the introduction of our earlier paper |J. The logarithm appearing in the 
time scale is a recent and interesting result of Bodeker ]/]], whose physical interpretation we 
discuss in ref. [BJ. 

Throughout this discussion, "hot" means that the temperature is large enough that the 
running coupling a(T) is small, that chemical potentials are ignorable, and that there is no 
spontaneous symmetry breaking. Examples of non-perturbative processes include chirality 
violation in hot QCD, and baryon number violation in hot electroweak theory (in its high- 
temperature symmetric phase). 

Bodeker @ has proposed an effective theory appropriate for the scales ( [LTD above. His 
effective theory is a classical field theory that involves only gauge fields with dynamics 
governed by the diffusive Langevin equation 

ffE = DxB-(. (1.2) 

Here, D is the covariant derivative acting in the adjoint representation. In Bodeker's pro- 
posal, £ is a Gaussian white noise random force, normalized asQ 

C(t,xKj(t',x')) = 2ag 2 T6 ab 6 l3 6(t-t')6^( X -x'), (1.3) 



where i,j and a, b are spatial vector and adjoint color indices, respectively. This effective 
theory is supposed to give a quantitative description of non-perturbative physics in the hot 
plasma to leading order in the logarithm of the coupling. In other words, corrections to this 
effective theory are suppressed only by powers of l/ln(l/g). In ref. ||, we showed that a 
can be interpreted as the color conductivity^ of the plasma, which is given by 

a (1.4) 

7g 



2 We will scale our gauge fields by a factor of g, so covariant derivatives contain no explicit couplings while 
the action (or energy) has an overall factor of 1/g 2 - In addition, we will take the gauge group generators, 
and also the gauge field A = AT, to be anti-Hcrmitian. Hence, the covariant derivative is simply 
D = V + A, and for the adjoint representation (T a )t, c — fbac with the structure constants f a bc real and 
totally anti-symmetric. 

3 We are using "color" as a descriptive name for some non-Abelian gauge field. It should be emphasized 
that all discussion of "color" is applicable to the dynamics of, in particular, the SU(2) electroweak gauge 
field. 
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where m p i is the plasma frequency and 



7g ^aC A T \n(l/g) (1.5) 

is the damping rate for hard thermal gauge bosons The ~ sign indicates equality 

at leading logarithmic order. [That is, we are not distinguishing ln(2/g) from ln(l/g) in 
Eq. |1.5| , but the coefficient of the logarithm is correct.] The plasma frequency m p \ is well 
knownQ at leading order in coupling and is of order gT. 

Bodeker's effective theory is well suited to numerical simulation because it is classical, 
insensitive to ultraviolet cut-offs 0, and when cast into A = gauge generates a straight- 
forward local equation of motion for the evolution of A: 

a ^-A = -D x D x A + C (1.6) 

A numerical investigation of Bodeker's effective theory and its implication for electroweak 
baryon number violation has been recently carried out by Moore |J. 

Nevertheless, there is something peculiar about the effective theory ( |1.2| ). In a high 
temperature plasma, static electric fields (or, more generally, longitudinal fields) are Debye 
screened [PJ. The screening distance is of order 1/gT, which is small compared to the 
spatial scale R ~ 1/ 1 g 2 T of interest to us. More generally, the longitudinal modes of the 
gauge field are screened while, at low frequencies to -C A;, the transverse modes are not. 
Because of Debye screening, it is the transverse electric and magnetic fields which are relevant 
for producing non-perturbative fluctuations at the scales ( [LTD quoted earlier. Longitudinal 
fields are irrelevant. Nonetheless, Bodeker's effective theory (|1.2|) does describe long-distance 
longitudinal fluctuations. The longitudinal fields are those pieces of E which contribute to 
D • E and which perturbatively correspond to polarizations parallel to the momentum k. 
Dotting D into both sides of (|1.2|) , one sees that 

crD E = -D C (1-7) 

D ■ E is therefore not zero, and so the fields in Bodeker's effective theory (|1.2|) are not purely 
transverse. Our purpose in this paper will be to show two things: first, that Bodeker's 
original derivation of the long-distance longitudinal dynamics relied on a questionable ap- 
proximation which ignored subtleties associated with longitudinal dynamics, but that the 
end result is nevertheless correct; and, secondly, that the longitudinal dynamics is irrel- 
evant and may be removed altogether, if one is purely interested in describing physical, 
gauge-invariant quantities that depend only on the transverse fields (for example, the rate 
of anomalous charge violation). The last point would be trivial in an Abelian theory, because 
then ( |1. 2D would be linear in the fields and could be projected into one equation involving 
only the transverse fields and another independent equation involving only the longitudinal 
fields. The point is much less trivial in a non-Abelian theory because of the non-linearity of 



4 For hot electroweak theory with a single Higgs doublet, for instance, = jg(5 + 2nf) g 2 T 2 at leading 
order in g, where m — 3 is the number of fermion families, and the adjoint Casimir Ca = 2 in ([L^). For 
QCD with n flavors of quarks, = 3 (l + §) 9 2 T 2 and Ca = 3, where n is the number of relevant quark 
flavors (u, d, s, c, b, t). 
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II. REVIEW OF BODEKER'S DERIVATION 



A. The effective Boltzmann-Vlasov equation 



We will refer to gauge fields associated with the scales of interest ( |1.1|) as "soft" fields. 



In contrast, the dominant excitations in the hot plasma correspond to momenta of order T 



and will be called "hard." On his way to deriving the effective theory (1^2) for the soft gauge 
fields, Bodeker pj first derived an effective Boltzmann-Vlasov equation for the interaction 
of those fields with hard excitations: 



[D t + v • D x ) W - E • v - £ = -5C[W], 



[2.1a) 



D U F^ U = J^ = ml (/lf(v)) v , (2.1b) 

where 

5C[W](v) = 7g (X(v, v') W(v')) v , , (2.2a) 

and 

T(v, v') = ^(v-v') - - (V ; x9 . (2.2b) 

V yjl ~ (V • V') 2 

Here, = Sm^ is the leading-order Debye mass (squared). The first equation fl2.1a|) is a 
linearized Boltzmann equation for the hard particles in the presence of a soft electromagnetic 
field, where W(v, x, t) represents the color distribution of those particles and v is a unit 
vector representing the hard particles' velocities.^ 5C represents a linearized collision term 
for 2— >2 scattering that randomizes the color charges of the hard particles |J, and £ is a 
source of random thermal noise. The second equation ( p.lbj) is Maxwell's equation, where 
all the fields on the left-hand side are to be understood as soft fields, and the current on 
the right side is the soft-momentum component of the current created by hard excitations. 
This current is proportional to the density W of hard particles and the velocities of those 
particles, where means (l,v). In the explicit form (|2.2j) for the collision term, (• • -) v 
denotes averaging over the direction of v and is a 5 function defined on the unit two- 
sphere with normalization 

5( 2 )(v-v')/(v')), =/(v). (2.3) 



See Refs. ]7||| for the derivation of the explicit form ( |2.2j b) of the linearized collision operator. 

One may avoid worrying about the details of noise terms such as £ until one reaches the 
final effective equation (|1.2|), at which point it is possible to then argue how the noise must in 



5 Technically, W is the adjoint representation piece of the density matrix describing the color charges of the 
hard excitations, summed over the various species of excitations and integrated over the energy of excitations 
(for a fixed direction of motion v) . It is normalized in a way that simplifies the resulting equation. See ref. 
J3[ for the explicit definition. 
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fact appear ||. However, since in this paper we will be discussing various subtleties, it will 
be useful to keep track of the noise explicitly at each step we consider. In particular, Bodeker 
|7]] derived that the appropriate normalization of the noise in the effective Boltzmann-Vlasov 
equation (|2.1| ) is related to the collision integral: 

(C(*,x,vK b (t',x',v')) = ^I(vy)5 ab 5(t-t')S^(^). (2.4) 

We will not review any further the origin of the effective Boltzmann-Vlasov equations ( [2 . 1| ) 
and direct the reader instead to Bodeker's original work and our alternative derivation 
|J. It is in the step from these kinetic equations to Bodeker's final effective theory ( |1.2| ) 
that subtleties in the treatment of longitudinal physics creep in, and that is the focus of this 
paper. 



B. Solving for W 



Bodeker obtains his final effective theory ( |1.2|) , at leading-log order, from the effective 



Boltzmann-Vlasov equations (|2.1|) by arguing that the covariant derivative terms in the 
Boltzmann equation ( p.laf ) are together of order g 2 T W and so can be ignored compared to 



the collision term, which is of order 'jgW ~ {g 2 T In g" 1 ) W and hence larger by a logarithm. 
There is an important subtlety to this approximation which will be examined in the next 
section. But accepting this argument at face value for now, if one drops the covariant 
derivative terms then the Boltzmann equation becomes simply 

E-v + f «<fC[W]. (2.5) 

Formally, the solution is 

W = (5Cy 1 (E-v + 0, (2.6) 

where 5C is to be understood here as an operator acting on the space of (adjoint- 
representation) functions of a unit vector v. This result for W yields the spatial current 
appearing in ( |2.1b| ), 



J = m^(v( ( 5C)- 1 (E-v + 0) v . (2.7) 

Next note that 5C preserves the parity (in v) of functions it acts on. In other words, 
SC maps even (odd) functions of v into even (odd) functions of v. (In contrast, the v • D x 
operator that we dropped does not.) Moreover, in the space of odd functions of v, 5C as 
given by (|2.2j ) reduces to simply 5C = 7 g . So, since 5C is a symmetric operator, and since it 
acts to the left on the odd function v in (|2.7|), we can replace (5C)~ l by 7 g 1 in that equation 
to obtain 

J = I^<v(E.v + 0>, = *E + C, (2.8) 

7g 

where 
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C = 3a<v£) v . (2.9) 

Using the correlation (|2.4j ) for £, one produces the correlation (|1.3| ) asserted earlier for £. 
Taking the spatial part of the Maxwell equations ( |2.1b| ) and dropping the dE/dt term which, 
for the scales (|1 . 1| ) of interest, is smaller (by four powers of coupling) than the crE term, 
one obtains Bodeker's final effective theory ( |1.2| ). 



III. LONGITUDINAL SUBTLETIES 

In the introduction, we noted that Bodeker's effective theory (|1.2p contains a fluctuat- 
ing longitudinal electric field. This may seem puzzling since longitudinal electric fields are 
Debye screened. In this section, we will take a closer look at how both Debye screening, 
and Bodeker's effective Langevin equation, do emerge from the Boltzmann-Vlasov equa- 
tions (pp . 



A. Zero mode of 5C 

In the last section, the Boltzmann equation for W was simplified, at leading-log order, by 
arguing that 5C dominates over the convective derivative D t + v ■ D x by a power of ln(g _1 ). 
This is not quite correct, however, because the operator 5C has an eigenvalue which is not 
order 7 g and which does not dominate over the convective derivative; specifically, SC has a 
zero mode. 

The necessity of this zero mode was noted by Bodeker, who observed that the effective 
Maxwell equation (|2.1b|) for the soft fields requires conservation D^J^ = of the current 



jfi = m'^ ) (v fJ- W( i v)) v generated by the hard particles. From ( |2.1a| ), this conservation requires 
(5C[W] ) v = 0, which is indeed satisfied by 



The fact (5C[1U]) V = can be rephrased to say that the symmetric operator 5C has 
null states: it annihilates anything that is independent of v. (This can be written in bra- 
ket notation in v-space as (constant |5C|W) = (lU|<5C|constant) = for any W .) This 
point will be important later on, so let us give an alternative way of understanding it. The 
collision term SC does not care, at leading-log order, about the dynamics of the soft fields. 
In particular, it does not care that the soft effective theory is a gauge theory, with a local 
color symmetry, instead of a non-gauge theory, with merely a global color symmetry. So, 
from the point of view of the calculation of SC at leading-log order, the theory could have 
been one where it was meaningful to talk about the total color charge of the system. If one 
then imagined adding an infinitesimal chemical potential u for this total color charge, the 
resulting equilibrium density would be 



n 



A-^Tl ~ = n + n (l±n )(3 gf x a T a + O(fi 2 ) (3.1) 



for each particle type, where no is the u = equilibrium distribution. In equilibrium, the 
collision term in a Boltzmann equation always vanishes by detailed balance. Different values 
of \i correspond to different equilibrium states, and the collision term must therefore vanish 
for all u. That means that the linearized deviation 



6 



Sn = n (l±n ) (3gn a T a 



(3.2) 



of the equilibrium distribution (|3.1| ) away from no must correspond to a null state of the 
linearized collision operator SC. The deviation ( |3.2j ) is isotropic and homogeneous — it is 
independent of both v and x. As a result, when re-expressed in terms of the function ^(x, v) 
used to parametrize color distributions of the hard particles in the linearized Boltzmann 
equation (|2.1a|) , the deviation ( |3.2|) corresponds to W^(x, v) = constant. That means that a 
constant W is a null vector of SC. But since collisions are local in x (in the effective theory), 
the x dependence of W is irrelevant, and so any W which does not depend on v is a null 
vector of the linearized collision operator SC. 



B. Longitudinal and transverse projections 

Before continuing, it is worthwhile to introduce longitudinal and transverse projection 
operators. Perturbatively, the longitudinal projection operator for the electric field is 

(ppert)y = fry = yi 1 yJ _ ^ 

V 2 

The gauge- covariant non-perturbative generalization is 

Pg^D^D*, (3.4) 

where D 2 means D • D. The transverse projection operator is of course 

p« = 5 ij - P[ j . (3.5) 

It is the longitudinal electric field which couples to external charges, since Gauss' Law reads 
D ■ E = p and since D • (Pp E) = 0. And it's the transverse electric field that is produced by 
D x B in the effective theory (|1.2|) since Pl (D x B) = 0. As mentioned in the introduction, 
the precise separation between longitudinal and transverse dynamics is not transparent from 
this simple discussion because of the non-linear dependence of D x B on the underlying vector 
potential A. 



C. Solving for W (again) 

To examine the difficulties caused by the presence of a zero mode of SC, we now return 
to the effective Boltzmann- Vlasov equation fl2.1| ) and will repeat the analysis of section [II B 



this time being more careful about how we treat the convective derivative compared to SC. 
Formally, the solution for W is 

t " = A + v 1 D. + fC (V - E + - (3 ' 6) 
or W = G(v ■ E + £), where G denotes the inverse of the linearized kinetic operator, 

G= [A + v • D x + • (3.7) 
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This fluctuation in the distribution of hard excitations produces a current response [from 
Eq. ( glbl) l of 



J = m 2 D (vG(v.E + ()) v , 
and the (color) charge density 

D-E = J° = m 2 D (G(vE + 0) v . 
One may easily check that the current is conserved (as it must be), since 

D J° + D • J = ml ((A + v • D) G (v • E + £)) 

= ml([l-5CG] (vE + O) 
= 0. 



(3.8) 



(3.9) 



(3.10) 



The (v • E + £) v term vanishes due to isotropy, (v) v = 0, and the lack of bias in the 
noise, (£) v = 0. And (SC G(v • E + £)) v vanishes because SC is acting (to the left) on its 
v-independent zero-mo de.[] 

D. The problem with the naive derivation 

In following subsections, we will discuss how to evaluate the operator inverse that defines 
G. To begin, however, it is useful to see how the zero-mode problem manifests itself in a 
simple example of matrix inversion. To this end, let us for the moment replace the Green 
function G of ( |3.7| ) by that of a simplified finite dimensional example. First, imagine that 
the gauge interactions are Abelian, so that D t and D x can be replaced by simply — iuj and 
ik, respectively. Next, imagine that the infinite-dimensional space of possible functions of 
v, on which 5C acts, is truncated to the four-dimensional space of functions that are either 
independent of v or linear in v. We wish to examine the matrix representing the action 
of —iuj + v ■ D x + 5C within this truncated space. In order to distinguish clearly between 
longitudinal and transverse physics, it is convenient to choose a basis {/ Q (v)}, a = 0, ... ,3, 
where 



/o(v) 



/i(v) = y/3ii ■ v 



(3.11) 



and e, are three mutually orthonormal unit vectors with e\ = k pointing in the direction of 



k. The overall normalization has been chosen so that (fifj) v 
elements of (fi\—iu + v ■ D x + 5C\fj) v are 



5jj. In this basis, the matrix 



7g 



\ 



-«*>+7g / 



(3.12) 



6 In slightly more explicit notation, this term is 7 g (X(v", v') G(v', v) [v • E + £(x, v)]) v , „. It vanishes 
because (X(v",v')) v „ = (Z(v", v')) v , =0. 
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The inverse operator, corresponding to G in our truncated space, is 



G 



trunc 



^k/(-i lg uj+lk 2 ) -iuj/(-i lg uj+lk 2 ) 



-ioj+j„ 



V 



\ 



;-z^+ 7g )- 1 j 



. (3.13) 



The uj —>■ limit is particularly simple: 



G 



trunc 



/ 37g 
' fc 2 

k 



V 



n/3 
k 





-1 



(3.14) 



7 S V 



In contrast, the naive derivation of Bodeker's theory corresponds to replacing — iu+iv-k+SC 
by SC, and the corresponding "inverse" would then be 



/ oo 



V 



7 E " 



(3.15) 



As one can see, there is no difference in the transverse sector (spanned by /2,3), but there 
is a huge difference in the longitudinal sector. As a particular example, consider the noiseless 
part of the current J, given in Eq. |3.8| ; namely mp(vGv • E) v . In the naive derivation, as 
represented by (|3.15|) , this contribution gives 



^d(v7 s 1 v-E) v 



3 7g 



E. 



(3.16) 



In our Abelian truncated-space calculation of G, however, it is instead given by 



ni 



^(G 11 P L + G 22 P T )E 



3 7g 



P T E. 



(3.17) 



in the u — > limit. The longitudinal part of E is projected out! This is quite different from 
the result of the naive derivation. 



E. Low-frequency, long-wavelength dynamics 

We will now show that, despite the major difference in G" 1 , one nonetheless does recover 
Bodeker's effective theory even for longitudinal physics. To do so, we will also return to the 
full, original non-Abelian problem and dispense with the truncated Abelian model of the 
previous section. 

If we restrict our attention to frequencies and wavenumbers which are small compared to 
the damping rate, u,k 7 g , then in the Greens' function G we could drop the convective 
derivative compared to 5C were it not for the fact that SC has a zero-mode. To deal with 
this, let Pq denote the projection operator onto the zero-mode of SC, so that 
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Po (/(v)) = </(v)> v , 



(3.18) 



and separate the convective derivative into zero-mode and non-zero-mode pieces, 

A + v • D = (A + v • D) P + P v • D + (1-Po) (A + v • D) (1-Po) • (3.19) 

To see this, note that A commutes with P , and that Pov • DP = 0. The last term of 
( |3.19j ) only perturbs the non-zero eigenvalues of 8C, and may be neglected provided uj and 
k are small compared to 7g . [For k = 0(g 2 T) this is a leading-log approximation.] The first 
two terms of ( p. 19 ) are rank one perturbations which will lift the zero-mode of the linearized 
kinetic operator. And one may evaluate explicitly the change in the inverse of an operator 
produced by adding a finite rank perturbation. In this case, a short exercise shows that 



G ~ [(A + v ■ D)P + P v ■ D + 6C]- 1 
= (1-Po) 5C-\1-P ) + 7 - 1 ( 7g - v D) 



Pn 



ry n. _ In 2 



( 7g -vD). 



(3.20) 



(To verify the last equality, recall that 5C is nothing but multiplication by 7g when acting 
on odd functions of v. Hence, 5C v • D P = 7g v • D P .) 

Now pause to note the correspondences of this result with the truncated Abelian results 
of the previous section. The full inverse ( |3.20| ) gives 



(G) v 
(vG) v 
(Gv) v 

(vGv) 



-7g 



7gA- 

7g A 



1 D 2 
3^ 



7-A - ±D= 



|D 2 

-1 



-1 



v ^ 7g < 1 



P 2 



1 

Dt^O 3 7g 



-Pt 



(3.21a) 
(3.21b) 
(3.21c) 
(3.21d) 



for uj, k 7g . These are simple non- Abelian generalizations of ( |3.13| ). (The factor of 3 
differences just reflect the normalizations of our chosen basis functions in the last section.) 

The form (|3.20| ) for G may now be inserted into Eqs. ( |3.8|) and (|3.9|) . For the charge 
density, one finds 



D • E = J° 



"7g 



7g A-|D 2 l (aD-E + D-C), 



(3.22) 



or 



A + (-D 2 + 



DE = -D-C 



(3.23) 



where a = |m5/ 7g , and C = 3cr(v^) v . And for the current, from (|3.8|) and (|3.22|) , 



J = aE + C 



mf 



D (D E) 



(3.24) 



Inserting this into the Maxwell equation —D t E + D x B = J gives 
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D t E + o-E 



a 



mf 



D(D E) = D x B-C 



(3.25) 



This local equation of motion is the exact result which follows from approximating G as 
shown in (|3.20|) . However, that form for G was based on the assumption that the fre- 
quencies and wavevectors of interest are small compared to the damping rate 7 g . Since 7 g 
is 0(g 2 Tln g^ 1 ), this means that uj is tiny compared to the 0(T/ Ing" 1 ) conductivity, and 
that k is much smaller than the O(gT) Debye mass. Hence, there is no point in retaining the 
D t E or D(D ■ E) terms in the effective equation ( |3.25| ). Dropping these terms immediately 
yields Bodeker's equation (1^). In other words, a more careful treatment of the effect of the 



zero mode in SC does not produce any difference (in leading-logarithmic approximation) to 
the resulting effective theory. 



F. Recovering Debye screening 

To see how Debye screening emerges from the kinetic theory fl2.1|) , return to Eq. ( |3.6| ) 
for W and now assume that the scales of interest are in the perturbative regime where 
k ^> g 2 T and/or u ^> g A T In g^ 1 . [In other words, the necessary conditions ( |l.lb| ) for 
non-perturbative fluctuations are not both satisfied.] In this regime, the gauge fields in 
the covariant derivatives appearing in the Greens' function G may be treated as smalLQ 
Expanding G in powers of the gauge field, the leading term, 



G~ [d t + y- V X + 5C] 



-i 



is diagonal in momentum space. Fourier transforming Eq. 
yields (to leading order in the gauge field) 

i k ■ E = m 2 D (G v) v ■ E + m 2 D (G£)v 

k E 

= m 2 D (G v • k) v — — + ml (G£)v 



(3.26) 

for the charge density then 



or 



with 



- k2 

-iml (1 + iuj(G) v 



k 2 + m 2 D (l+ iw<G) v )] i k • E = m 2 D k 2 (G^v 



G(u, k) = [ — iuo + iv ■ k + 5CY 1 



(3.27) 



(3.28) 



(3.29) 



and E and £ now denoting the (u, k) Fourier components of these fields. In the first step of 
( [3.271 ), we used the fact that, with gauge fields neglected in G, the only vector which (Gv) v 



7 For a more detailed justification, based on a computation of the power spectrum of gauge field fluctuations, 
see Rcf. ||. ' 
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can depend upon is k, and therefore only the longitudinal component of E can contribute 
to the result. The following step used (G (zv • k — iuj)) = 1 — (G8C) = 1 which is another 
consequence of the zero mode in 5C. The result ( |3.28| ) shows that k • E satisfies a diffusive 
Langevin equation in which the noise and damping depend on (G£) v and {G) v , respectively. 
The power spectrum of charge density (or D • E) fluctuations is defined as 



Pl(^? k) = dt <i 3 x e 



Auit— ik-x 



(J°(*,x)J°(0,(%. 



(3.30) 



Using (|3.28|) to express D • E in terms of the noise £, and then recalling that the covariance 
of the noise £, as given by ( |2.4| ) and (2^), is proportional to 5C, allows one to write the 
power spectrum as 

p L (u, k) = m 4 k 4 ((G'e)v(^t} v ,) £ /\k 2 + m 2 D (l + iu(G) 

= 2g 2 T ml k 4 Re (G) v / | k 2 + ml (l + iu(G) 
g 2 k* 



2T 
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Im 



k 2 + m* (l + zu;(G% 



(3.31) 



Once again, this answer is valid provided k ^> g 2 T and/or uj ^> g 4 T\n g" 1 since we neglected 
(for this discussion only) the soft gauge fields in the covariant derivatives. Furthermore, k 
and io must be small compared to T, since this is a basic requirement for any kinetic theory 
description to be valid. 

Linear response theory, applied to the underlying quantum field theory, shows that the 
power spectrum 



PL(co>,k) = / dtd 3 xe 



Awt— ik-x 



i{j°(t,x),J°(0,0)} 



(3.32) 



(defined with a symmetrized ordering of quantum operators) is related to the retarded 
charge- density charge-density correlator 



D R (u,k) 



dt g? 3 x e 



iujt— ik-x 



0(t)( J°(*,x),J°(0,0) 



by the fluctuation-dissipation relation 

p h (u,k) = [2n(u) + 1] lmD R (oo,k) 



(3.33) 



(3.34) 



where n(u) = \er w — l]^ 1 is the equilibrium Bose distribution function. And the retarded cor- 
relator -Dr(u;, k) is the analytic continuation of the Euclidean space (time-ordered imaginary- 
time) correlator D E (iu n ,k) from the imaginary Matsubara frequencies to (just above) the 
real frequency axis. 

By comparing ( |3.31|) to the form ( |3.33|) , the kinetic theory result for the power spectrum 
may be converted to an equivalent result for the retarded correlator. The leading factor of 
T juj is just the low-frequency (classical) limit of the Bose distribution function, and hence 
the retarded charge density correlator is 



S 2 k 4 



D R (u;,k) = g 2 k 2 



k^ + ml (l + iw(G) Y ) ' 



(3.35) 



12 



The local (and temperature-independent) g 2 k 2 term, which does not contribute to the imag- 
inary part, is determined by the current-current Ward identities, or equivalently by the 
requirement that remain bounded as k — > oo. The low frequency limit, 



Al(0,k) 



g 2 k 2 m 2 } 
k 2 + rap 



m 



m D - 



D 



k 2 + mn 



(3.36) 



reproduces the correct static equilibrium Debye-screened charge density correlations.^ More 
generally, the kinetic theory answer ( |3.35| ) reduces to the known hard-thermal-loop result 
whenever the frequency or momentum is large compared to the damping rate j g . In this 
domain, the details of the scattering operator 5C are irrelevant, and G(o;,k) may be ap- 
proximated by i[u — v • k + ie]^ 1 . The resulting average over v may then be performed 
analytically, and yields the standard HTL result for the self-energy. 

If neither k nor uj are large compared to the damping rate 7 g , then the detailed form of 
SC is significant. Evaluating G(u, k) is non-trivial if k is comparable to the damping rate. 
However, if k and u are both small compared to 7 g , then the previous representation (|3.21a ) 
may be used. Perturbatively, it gives 



<GV,k)>, 



7g 



±k 2 - ZCJ 7g ' 



(3.37) 



which gives 



PiXw,k) 



2g 2 T 



k 2 . 



a 



(3.38) 



This is the same result which emerges directly from the relation o\DE = — D-£ in Bodeker's 
effective theory, combined with (|1.3j ) for the noise variance. This result is valid in the overlap 
of the perturbative and damping-dominated domains, that is when \k 2 — iua\ 3> g 4 T 2 and 



IV. IRRELEVANCY OF LONGITUDINAL DYNAMICS 

We have seen that Bodeker's effective theory 

(tE = DxB-C, (4.1) 

within its domain of validity, does correctly describe both longitudinal and transverse fluc- 
tuations. However, it is the transverse part of the gauge fields which are responsible for in- 
teresting non-perturbative phenomena such as topological transitions and associated baryon 
non- conservation. One may wonder if it is possible to formulate an equally valid effective 



8 The factor of <? 2 is present merely because we have chosen to scale all our gauge fields by g relative to 
the usual perturbative conventions. In, for example, Coulomb gauge, Z?r equals the one-loop gauge field 
self-energy IT 00 (which is just nv^ in the static limit), plus the one-particle reducible contributions which 
sum to n 00 D 00 n 00 , where D oa = (A A ) = -l/(k 2 + II 00 ) is the Debye-screened A propagator. 
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theory which describes only transverse physics. On the face of things, this should be easy; 
just insert a transverse projection operator to eliminate the longitudinal part of the noise, 



aE = D x B - P T C • (4.2) 

This produces an effective theory with no longitudinal dynamics whatsoever, PlE = 0. In 
the case of an Abelian theory, the trivial decoupling of transverse and longitudinal parts 
of the gauge field would make it obvious that Bodeker's theory fl4.1|) and the transverse- 
projected theory Q4.2| ) describe exactly the same transverse dynamics. But for our non- 
Abelian theory, the dependence of covariant derivatives and projection operators on the 
gauge field makes this decoupling far less obvious. The remainder of this paper is devoted 
to showing that it is almost true that equations ( |4.1| ) and ( |4.2| ) do in fact generate iden- 



tical transverse dynamics. The "almost" caveat reflects subtleties associated with the fact 
that white noise cannot be considered smooth even on infinitesimal time scales. Investigat- 
ing these subtleties will reveal that the naive transverse equation ( f4.2|) must be corrected 
but can then be made exactly equivalent to Bodeker's equation ( |4.1| ) with regard to the 
transverse dynamics. It should be emphasized that this investigation was motivated by 
theoretical curiosity, not practical convenience .0 Bodeker's equation ( [lip is local, whereas 
the transverse-projected equation ( |4.2j ) is not. Consequently, numerical simulations are far 
more easily performed in the original theory than in any transverse-projected variant. 

To be precise, by "transverse dynamics" we refer to all physical observables that do not 
depend on A when expressed as gauge- invariant functions of the fields A^. The magnetic 
field B = D x A does not, of course, depend on A . Writing E = DA — (dA/dt), it is 
easy to see that the transverse electric field PpE = — Pt (dA/dt) does not either. Conse- 
quently, an example of a physical quantity which depends only on the transverse dynamics 
is the topological charge (or change in Chern-Simons number) of the gauge field, which is 
proportionaQ to /tr[E • B] = / tr[(P T E) • B].Q 



A. Naive equivalence 

We first wish to paint with a broad brush. We will for the moment ignore all subtleties 
and discuss how, if one implicitly and incorrectly (and only when advantageous) treats the 
noise C(x, t) as a smooth function of t, one may show that the two theories ( f4.1|) and ( |4.2| ) 
should generate the same transverse dynamics. We will wait until section [IV B| and its sequel 



9 In fact, historically, it took us some time to recognize that longitudinal physics for u> <C k <C 7 g is 
correctly reproduced by Bodeker's effective theory, and the discussion in this section was motivated by the 
desire to show that it simply doesn't matter. 

10 The precise formula is AN(t) = (1/8tt 2 ) J* dt d 3 x EfBf. 

n The topological transition rate (or Chern-Simons number diffusion constant), is an important ingredient 
in scenarios of electroweak baryogenesis. Understanding the applicability of numerical simulations using 
Bodeker's effective theory for extracting the topological transition rate motivated this investigation. See |I 
for such recent numerical work and related discussion. 
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to correct this discussion by taking into account the non-smooth nature of Gaussian white 
noise. 

It is simplest to initially consider both theories ( |4.2| ) and ( |4.1| ) in A Q = gauge: 



a— A = — D xB + ( 

at 



(4.3) 



and 



a— A = -D x B + P T C • 

at 



(4.4) 



For the moment, imagine a particular instantiation of the white noise C( x )^) — that is, 
consider a particular member of the Gaussian ensemble of noise functions. Suppose that 
Bodeker's equation ( f4.3|) is satisfied by a gauge field A(x, t). Now rewrite Bodeker's equation 
in the form 



a(j t A-a- l P^C] =-DxB + P T C 



(4.5) 



Using the explicit form ( |3.4|) of the longitudinal projection operator, this can be written as 



<r(|A-D4 



-DxBl P T C 



(4.6) 



where Aq is simply a (suggestive) name for 



(4.7) 



For a particular instantiation of the noise (and any initial condition), the solution to 
JD may be interpreted in two different ways. On the one hand, A is, by construction, an 
A = gauge solution to Bodeker's equation (|4.1|) . On the other hand, if one says that A is 
actually the time component of the gauge field, then the left-hand-side of ( |4.6|) is just — crE. 
Therefore, A^ = (A , A) is a solution to the projected equation ( |Q| ) in the particular gauge 
where A = a^D^D ■ £. 

But, given a solution (A , A) to ( f4.2|) with A ^ 0, one may always gauge transform 
back to A = gauge. The result will be a gauge field A which obeys 



d X 
a— A 

dt 



-DxB + ?tC 



(4.8) 



This is just the A = transverse equation (|4.4j), except that the noise has been gauge 
transformed by the transformation which takes (Aq, A) into Aq = gauge: 



<f = u ab c\ 



with 



U = Texp 



A n dt 



Texp 



a 



* 1 
— B-Cdt 

o JJ Z 



(4.9) 



(4.10) 
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where T signifies time ordering, with the latest times on the right. 

The distinction between £ and its gauge transform £ will not matter, and our two theories 
( |4.3| ) and ( f4.4| ) will be equivalent (subject to earlier caveats), provided the distribution 



£ u = U ab C, is Gaussian white noise, just like the distribution of the original £. If our 
transformation U was not a function of the noise, this would be trivial because then we 
would have 



;C C ) = u °° u \CC ) • (4-11) 
Since the £ correlator is proportional to 5 cd , and since u ac jj hd § cd = § ab ^ it would follow that 

c a c7 = (cc7- ( 4 - 12 ) 



Even though our transformation U depends on the Gaussian white noise £, this result naively 
remains true. Consider, for instance, the equal time correlation 

C a (t)C\t)) = (u a \t)U bd {t) C(t)C d (t)) . (4.13) 



Because the noise correlation is local in time, while U (formally) depends only on the noise 
prior to t, this can be factorized:Q 

C a (t)C\t)) = (u ac {t)U bd {t)) (C(t)C d (t)) . (4.14) 

The C correlation is again proportional to 6 cd , which again contracts the £/'s and eliminates 
them, so that we arrive at Q4.12|) as desired. A similar argument shows that unequal time 
correlations of £ vanish, as they should. 



B. So what's wrong? 

1. A toy model 

To see what goes wrong with the previous equivalency argument, and to understand 
what it has to do with the short-time nature of white noise, it is instructive first to consider 
a system much simpler than non-Abelian gauge theory. Forget about field theory and 
instead imagine stochastic dynamics of a classical particle moving in two dimensions in a 
rotationally-symmetric potential V(r): 

j t r = -VV + C (4.15) 
(d{t)C j (lf)) = 2T6 ij 6(t-lf). (4.16) 



12 In fact, the dependence (or lack thereof) of U(t) on the noise at exactly time t is ill-defined and depends 
on the details of time discretization, as discussed in the next sub-section. 
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For convenience, we have normalized the analog of cr to 1. Imagine also we care only about 
the radial dynamics of this system and not at all about the angular dynamics. 

Comparing to the gauge theory problem, r above is analogous to A in A = gauge, the 
radial dynamics to the transverse gauge dynamics, and the angular dynamics to the longi- 
tudinal dynamics. Circles about the origin are analogous to gauge orbits of 3-dimensional 
gauge configurations under 3-dimensional gauge transformations (since, in the gauge theory, 
infinitesimal displacements in the longitudinal direction are of the form AA(x) = DA(x), 
which is the form of an infinitesimal 3-dimensional gauge transformation). Eq. ( |4.15| ) is 
analogous to Bodeker's effective theory Q4.3p, and the analog of the transverse-projected 
theory ( |4.4| ) is then 

^ r = _ W + P r C, (4.17) 

where the radial projection operator P r is 

P ij = = S ij _ Qi§3_ ( 41g ) 

Just as in section |1V A] , we can make a sloppy and not quite correct argument that the 
unprojected equation ( |4.15| ) and the projected equation ( |4.17|) are equivalent. A transforma- 



tion from a solution r of the unprojected equation to a solution r of the projected equation 
appears to be 



or equivalently 



Ur, C = UC, (4.19) 



U-'f, C = U~ l C, (4.20) 



where, if r and £ are represented by complex numbers r x + ir y and ( x + i( y , U can be written 
in a form quite analogous to Q4.10| ): 

U = exp(-i£§-Cdt\. (4.21) 

U simply rotates away the accumulated motion in the angular direction, so that the projected 
motion, at every instant in time, becomes purely radial. Naively plugging fl4.20| ) into the 
unprojected equation ( (4.15| ), and implicitly but incorrectly assuming that £ is a smooth 



function of time, yields the naive projected equation ( f4.17|) for r. 

One can immediately see that the two equations ( |4.15| ) and ( |4.17| ) cannot, however, 



actually describe the same radial dynamics. In the unprojected case ( 4.15 ), there is zero 
probability that the system would ever pass exactly through the origin r = 0. The projected 
case ( |4.17| ), however, just describes one-dimensional motion, parameterized by r, along a line 
of constant 9. That is, we could fix 9 and just replace (|4.17| ) by the one degree of freedom 
equation 
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(There does not appear to be an analog of this simplification in the gauge theory; see 
Appendix ||.) As long as there are no infinite potential barriers, this one-dimensional system 
will eventually fluctuate through any value of r, including r = 0. 

To understand the discrepancy, we need to properly understand the small time behavior 
of white noise Langevin equations such as (|4.15|) and ( |4.17| ). The standard way of defining 
what such equations actually mean is to discretize time and only at the end of the day take 
the continuous time limit. 



2. Time discretization ambiguities 

Before proceeding, we have to dispose of an instructive red herring concerning the time 
discretization of our various stochastic equations. It is not always true that continuum-time 
stochastic equations like the ones we have been writing down have an unambiguous meaning. 
To understand the possible ambiguities, imagine that instead of being interested in only the 
radial dynamics of our toy model, we were instead interested in only the angular dynamics, 
and so had proposed a projected equation of the form 

|r-a<, (4.23) 

Ptf = 0*0* = &> -fV J '. (4.24) 

This continuum equation appears to describe motion for which the radius r remains constant. 
Now imagine discretizing time with small time steps of size e, so that 

e- 1 Ar = P C, (4.25) 

(d(t) Q(t+ne)) = 2Tr l 5 i3 5 n0 . (4.26) 

There is an ambiguity in the schematic way we have written the discretized equation ( [4.251 ): 
we have not made it clear whether the direction 9 implicit in Pg is supposed to be evaluated 
at the starting point of the tiny time interval, the end point, or somewhere in between. In 
the first case, the value of r will drift out a little bit, as in fig. |l]a. In the second case, it will 
drift in a little bit, as in fig. [I]b. If we pick a symmetric convention, where we evaluate 9 at 
the midpoint, then r will remain constant, as in fig. |I]c. 

In non-stochastic equations, such discretization choices become irrelevant in the con- 
tinuum limit e — > (though they may have significance for the practicality of numerical 
calculations). For stochastic equations, however, the e — > limit is more subtle because, by 
( |4.26| ), the amplitude of the white noise £ is order 




(4.27) 

and diverges as e — > 0. The drift Ar in figs. |l]a and b is therefore of order 

± Ar ~ yr 2 + (eC) 2 - r ~ „ e 2L (4.28) 
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FIG. 1. Radial motion generated by ( 4.25 ) depending on whether 9 is defined by the (a) beginning, (b) 
end, or (c) middle of a discrete time jump t —> t + e. 



for a time interval e. That means that the drift per unit time, Ar/e, is finite as e — > 0, and 
so the continuum limit really depends on one's discretization conventions. 

In the unprojected equation (|4.15| ), there is no such discretization ambiguity. And in 
our actual toy model equation ( |4.17| ) with radial projection, there is no such ambiguity 
because motion in the f direction, unlike in the 9 direction, is straight — r does not change 
between one end of the interval and the other. The situation is slightly more complicated 
for the transverse-projected equation (|4.4j) for gauge theory, however. There, motion of A 
in the transverse direction does change the transverse projection operator Pj- However, we 
demonstrate in section [IV C 2| that this change turns out to be high enough order in 5 A that 
discretization ambiguities do not arise. 

The upshot is that the continuum stochastic equations ( |4.17|) and Q4.4J ) for the radial- 
projected toy model and the transverse-projected gauge theory are not ambiguous. However, 
we shall next see that the very same discretization issues affect the transformations we used 
to argue that they were equivalent with their unprojected counterparts. 



3. Centrifugal drift 



The way we proposed obtaining the projected equation ( [4.171 ) from the unprojected one 
( |4.15| ) was by rotating away the accumulated motion. Imagine a single time step of the 
discretized unprojected equation. Then0 



r(t + e) = r(t) - eW(r(t)) + 60). 
The motion of the radial coordinate r is then 



(4.29) 



13 Whether it's V(r(t)) or V(r(t+e)), or the average of the two, does not matter in the continuum time 
limit e — > 0. 
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r(t + e) = |r(t) - eW(r(t)) + eC(t)\ 



r-e1/'(r) + er-C + — Kl 2 

2r 



+ 0(e 3 / 2 ). 



(4.30) 



Given that, as e — > 0, a large number of successive tiny time steps will occur before the 
system appreciably changes position, the positive |£| 2 term can be replaced by its statistical 
average ( |4.26| ): 



r(t + e) 



r — eV'(r) + er - ( 



eT 



(t). 



(4.31) 



The distribution of r ■ £ does not care about the direction of f , and ( f4.31|) can be rewritten 

as 



where ( r is uncorrelated white noise and 

V<&{r) = y{f) — Tlnr. 



(4.32) 



(4.33) 



The continuum projected equations that are truly equivalent to the unprojected one are 
therefore ( 4.17]) or ( 4.22 ) with V replaced by V e g. The addition of the lnr term in ([4.33 ) 



now provides a "centrifugal potential" which prevents the one-dimensional system ( [4.22 
from passing through r = 0. 



4- Equilibrium distributions 

The exact form of the "centrifugal" correction was really determined from the very start. 
As we shall briefly review in section [IV C4| , the equilibrium distribution in r generated by the 



unprojected equation fl4.22|) is proportional to exp(— V/T). That means that the probability 



distribution for the radial variable r must be proportional to 

2vrr exp(-V/T) oc exp(-V cfr /T), (4.34) 

since the 2irr is the volume of the symmetry orbit. But exp(—V c g/T) is precisely the 
equilibrium distribution generated by the projected equation (|4.22|) if V is replaced by V e g. 



In the gauge theory case, there will be no analog to the one-dimensional radial equation 
( |4.22j ), and so it is worthwhile to understand how the equilibrium distribution could have 
been deduced directly from the two-dimensional projected equation (|4.17 ). This requires 



deriving the Fokker-Planck equation that is associated with a given Langevin equation and 
which describes the time evolution of probability distributions V(r). (This will be discussed 
explicitly in the gauge theory case below.) For the naive projected equation (|4.17| ), one finds 
that 

P(r) oc ^tXIIl (4.35) 
zrrr 

is a (two-dimensional) equilibrium solution. But if we correct the naive equation by replacing 
V Kff> then we indeed recover the unprojected result 
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exp(-V eS /T) 

V(r) oc = exp(— V T) (4.36) 

27rr 

as the equilibrium distribution in r.Q [Do not confuse the two-dimensional distribution 
( |4.36| ) for r with the one-dimensional radial distribution ( |4.34j ) for r. Both describe the 



same equilibrium ensemble.] 



C. Gauge Theory 

1. Time discretization ambiguities 



The transverse-projected version (f4.4| ) of the soft effective theory is a Langevin equation 
of the form 

j t Qi = -d t V(q) + e ia (q) (a, (4.37a) 

(Ut)( b (t')) = 2T6 ab 5(t-t'), (4.37b) 

where, for the moment, we are using notation natural for a system with a finite number 
of degrees of freedom. In the field theory case, the dynamical variables q are the values of 
the gauge fields at different points in space and <9j becomes a functional derivative 5/5 A. 
The functions ej a (q) characterize the coupling of the noise to the dynamical variables q; for 
the gauge theory this is the transverse projection operator Pt (which depends on the gauge 
field A). 

To define exactly what is meant by this equation, imagine discretizing time into very small 
time steps of size e.Q Stochastic equations of the form ( |4.37| ) are generically ambiguous if 
the coupling e ia (q) to the noise has non-trivial dependence of q, because of the ambiguity, 
discussed earlier, of when to evaluate q. In the discretized equation, 

«(' + ')- «M =HW , + efc (J w (4.38) 



(CWa«')> = 2Te- 1 ^fe, (4.39) 



14 Because the projected equation conserves 8, there is a family of two-dimension al eq uilibrium solutions to 
the projected toy model equation: namely, the rotationally invariant distribution ( |4.35| ) can be multiplied by 
an arbitrary angular distribution f{9). This non-uniqueness is, of course, irrelevant if one is only interested 
in rotationally invariant observables. The appearance of an arbitrary function of 9 in the general equilibrium 
distribution is a reflection of the non-ergodicity of the projected two-dimensional evolution equation. As 
discussed in Appendix 0, for the transverse-projected gauge theory there does not appear to be any analog of 
a conserved gauge-orbrc coordinate 9 and, so far as we know, the transverse-projected gauge theory remains 
ergodic. 

15 The following discussion roughly follows the presentation in sections 4.7 and 4.8 of ref. |l2), although 
our normalizations are slightly different. 
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this ambiguity appears as dependence on a parameter a which controls the time at which 
the right-hand side is evaluated. For example, a = corresponds to a forward time deriva- 
tive and is known as the Ito convention, a = ~ to the symmetric derivative, known as 
Stratonovich convention, and a = 1 to a backward time derivative. The precise meaning of 
evaluation at time t + ae is to expand in ae. Keeping in mind that the amplitude of the 
noise ( is e~ 1//2 , and using the equation of motion itself, the terms in the expansion which 
are non-negligible when e — > are 



qi(t + e) - qi(t) 



(4.40) 



-diV + e ia Ca + a e (dj 

The product ( a (b may be replaced by its expectation, giving the final discretized equation 
qi(t + e) - qi(t) 



-d i V + 2aT(d j 



(4.41) 



The term proportional to a is a convention-dependent "drift" force. The naive continuous- 



time formulation (4.37) does not, in general, uniquely specify the dynamics. 



2. Vanishing ambiguity 

We shall now show that the ambiguity vanishes for the transverse-only noise equation of 
( |4.4|) . This implies that there is no real issue of convention dependence for this application. 



We work in continuous space (rather than working on a spatial lattice, which would be 
more relevant to numerical simulations but also more complicated) . The degrees of freedom 
in the gauge theory case are labeled by coordinates X = (x.,i,a), where i is a vector index 
and a an adjoint color index. It will be convenient to introduce combined labels for several 
different choices of coordinates: 

X = (x,i,a), Y = (y,j,b), Z=(z,k,c). (4.42) 

The noise coupling e ia introduced above becomes the (matrix elements of the) transverse 
projection operator 

P XY = 5 l] 5 ab <J(x-y) - {p i D- 2 D j )%. (4.43) 

This operator is symmetric in X and Y, and the drift force discussed earlier is proportional 
to 

Pxzj^Pzy- (4.44) 
d A x 

When taking the variation of Pzy, the variation must act on the left-most covariant deriva- 
tive in ( 4.43Q , since otherwise that derivative will annihilate against the Pxz factor. One 



thus obtains 

Pxz t^- Pzy oc / P X z S lh T a ce <f(x - z) (D- 2 D^ b 

OAx Jz V /zy 



oc dT a ae (p-^jy^f - {D i D~ 2 D i )^TZ B (p- 2 D j ) e ^ (4.45) 



22 



in d spatial dimensions, where no integration over x is implied. The first term vanishes 
because the adjoint generators are anti-symmetric in (abc) and so T£ e = 0. The second 
term vanishes because (D l D~ 2 D l )'£ x is symmetric in (ac) and so vanishes when contracted 
with the anti-symmetric T£,. So the convention-dependent drift force is exactly zero 



3. Centrifugal drift 

We now return to the transformation, between Bodeker's effective theory and the 
transverse-projected theory, in order to derive the gauge theory analog of the centrifugal 
correction discussed for the toy model in section |IV B 3| . The time-discretized version of 
Bodeker's effective theory is 

A(t+e) = A(t) - e - (f^(t)-C(t)) , (4.46) 



where 



V=\g~ 2 B?B? = \g- 2 F%F° (4.47) 



is the potential energy associated with the magnetic field, the noise correlation is 

(C(t, x) $(t+ne, x')) = 2g 2 Tae- 1 5 ab 5 lJ 5 n0 5(x-x') , (4.48) 

and we have chosen to use a forward time difference. We now want to apply (a discrete 
version of) the gauge transformation U that was introduced in the time-continuum case, 
( |4.10|) , to eliminate longitudinal motion of A. For simplicity of presentation, we shall focus 
on one single, time step from t to t+e and discuss how to transform A(t+e) relative to 
A(t) in order to eliminate the longitudinal motion introduced during that step. Consider a 
gauge transformation U which equals the identity at time t, but is a non-trivial infinitesimal 
transformation at time t+e, 

U{t, x) = 1, U{t+e, x) = exp a(x) . (4.49) 

For the moment, we leave the generator of the transformation, a, arbitrary. The gauge 
transformed field is 



16 The drift force is also proportional to T which, in the field theory case, is T8(0). If we were only 
interested in perturbative physics, we could have chosen to work in dimensional reguiarization, which sets 
(5(0) to zero. 
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In the general case, a suf ficient condition for the ambiguity to vanish can be expressed as follows. Suppose 



the potential V(q) of ( 1.371 ) is invariant under some symmetry transformations that have the infinitesimal 
form q — > q + A0 a (q), where a indexes the independent symmetry generators (and A is infinitesimal). 

Define a metric g a/3 = 6 a ■ 8^ on the space of symmetry generators. As in general relativity, let g a p denote 
the (matrix) inverse of the metric. Now suppose that the noise coupling equals the projection operator 
e ij — % ~ Sfg a p6j. One can then show that the ambiguity (djei a )ej a vanishes if di6" is anti-symmetric 
under interchange of i and j. This anti-symmetry condition is indeed satisfied in both the radial-projected 
toy model and transverse-projected gauge theory. 
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A = U(V + A)U- 1 . 
Expanding in powers of the generator a at time t+e, this gives 
A(t+e) = A(t+e) - Da + ~[Da, a] + 0(a 3 



(4.50) 



(4.51) 



where Da = Va + [A(i+e),a]. Using the equation of motion (|4.46|) to rewrite A(t+e) in 
terms of A(t) gives 



A{t+e) = A{t) 



a UA 



C 



Da 



a 



5A 



C,a 



+ |[Da,a] + 0(a 3 , Tea 2 ), (4.52) 



where now all the covariant derivatives involve the gauge field at time t. Choosing the 
infinitesimal generator to equal 



a 



a 



ZT 2 D ■ C 



(4.53) 



as implied by our previous discussion [c.f. ( |4.10| )1, will cause the Da term to cancel the 
longitudinal projection of the noise £. Since the noise £ is of order e -1 / 2 , this means that 
a is of order y/e. We need to keep all terms in ( |4.52| ) which are of order e. The term 
(e/cr) [5V, a] is 0(e 3//2 ) and may be neglected. Consequently, 

A(t+e) = A(t) + i(-^L + P T A -4k" iDZT 2 D • <, D' 2 !) ■ C] + 0(e 3 / 2 ). (4.54) 



a V 8 A 



<j- 



Once again, we can replace the terms quadratic in noise by their statistical averages, as 
given by (|4.48 ). After some manipulation, one finds that this yields 



A a (t+e,x) = A a (t,x) + - 



5V 



Te 



P T C(t,x.)\+—f abc (DD 



a \ 5A a (x) 



1-2\cb 



a 



+ 0(e 3 



A a (t,x) + 



5V, 



(•ir 



+ P T C°(t,x) +0(e 3 



a \ <5A a (x) 
where the effective potential V e ff is 

V eff = V- iTTrln(-D 2 ) = V-Tln y / det(-D 2 ) . 



(4.55) 



(4.56) 



As shown in Appendix ydet(— D 2 ) is the volume of the gauge orbit containing a given 
spatial gauge field. Consequently, this logarithmic correction to the potential is completely 
analogous to the "centrifugal" potential appearing in the rotationally invariant toy model. 
The upshot is that the projected equation which is actually equivalent to Bodeker's effective 
theory differs from the naive projected equation (|4.4|) by the replacement of V by V e ff : 



(4.57) 



a— A = — D x B + IT— Tr ln(-D 2 ) + P T C 
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4- The equilibrium distribution 



It's interesting to examine what happens if one converts a Langevin equation of the 
generic form (|4.41| ) into a Fokker-Planck equation for the evolution of the probability dis- 
tribution P(q, t). One finds (see for example |12[) that 



^-V = d i [Td j (e ia e ja V) + {d.V - 2aT [d 3 e ia )e ja }v\ (4.58) 

If e ia (q) were just S ia , as in Bodeker's equation (fOf), then the equilibrium distribution (the 
solution to dV/dt = 0) would simply be V = exp(— V/T) up to an overall normalization 
constant. 

For the naive radial-projected toy model Q4.17D with Ci a = fif a , 



dj(e ia e ja ) = — = f~ ln(2vrr), (4.59) 
and (djei a )ej a = 0. This leads to the equilibrium distributions (|4.35| ) quoted in section 

[vbj . 



In the gauge theory case, we have seen that the convention-dependent drift force vanishes 
for the transverse-only noise equation (|4.4|) , but there still remains ei a dependence in the 
Fokker-Planck equation. Plugging in the transverse projection operator for e, one finds 

S S 
dj(e ia e ja ) ^—P XY = \P XY trln^D 2 ), (4.60) 

OAy OAy 

and, solving for the equilibrium distribution, 

/det{-D 2 ) 



up to an overall normalization constant.^ As mentioned above, ydet(— D 2 ) just represents 
the gauge orbit volume, and the above distribution is analogous to the toy model result 
( |4.35| ). As with the toy model, however, if we examine the transverse theory that is truly 
equivalent to the unprojected theory, then we should replace V — > V e s and we recover the 
correct equilibrium distribution 

V = eXp( ~ Vcff/T) = exp(-V/T) . (4.62) 
/detf-D 2 ) 



18 In the more general notation of footnote [l7[ the assumption that di6" is anti-symmetric in i and j leads 
to dj(ei a ej a ) — ^Sijdj In yj~g and V = exp(— V/T)/ ^/g, where g is the determinant of the inverse metric g a /3- 
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APPENDIX A: THE VOLUME OF GAUGE ORBITS 



The natural metric on the space of gauge field vector potentials is 



ds 2 



tr 



dA f • dA 



(Al) 



This is the unique metric (up to an overall multiplicative constant) which is invariant under 
both gauge transformations and spacetime symmetries. The gauge orbit passing through a 
particular gauge configuration A consists of all gauge transforms of A. Within a neighbor- 
hood of A, configurations on the gauge orbit may be parameterized as 



e- A De A . 



(A2) 



where A is an arbitrary generator of the gauge group (i.e., A(x) = A a (x) T a is a Lie-algebra 
valued function of x), and D = V + A is the covariant derivative with gauge field A. Since 
A A — A ~ DA, the induced metric on the gauge orbit, evaluated at A, is just 



ds 2 



orbit 



J tr [(D5A) t • (DM)] = j tr [5A\-D 2 ) 5A 



(A3) 



Consequently the induced volume element on the orbit, evaluated at A, is 



dv = Jdet(-D 2 ) dA, 



(A4) 



where dA = n x ,a^°(x) denotes the flat measure on the gauge algebra. But the gauge- 
invariant Haar measure on the gauge group is just dA when evaluated at the identity. And, 
because the functional determinant det(— D 2 ) is gauge invariant, it is constant over the 



gauge orbit. So, globally, the volume element dv equals ydet(— D 2 ) times the Haar measure 
on the gauge group. Hence, 



orbit volume 
gauge group volume 



det(-D 2 ) 



1/2 



(A5) 



and so y / det(— D 2 ) is the gauge orbit volume up to an overall A-independent normalization 
factor. 
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APPENDIX B: NO GAUGE THEORY ANALOG TO TOY MODEL 9 



Return, for a moment, to the toy model described in Sec. [IV B| . The projected equation 



( [4.171 ) h as two obvious properties. First, the particle always moves in a direction perpendic- 



ular to the gauge orbits r = const. Second, moving according to this equation, the particle 
cannot reach any point in the configuration space, but instead remains confined to a slice in 
a configuration space where 9 is a constant. In particular, starting from a point (r, 9), one 
cannot reach any point that is gauge equivalent to it, except the original point. In other 
words, if one fixes the gauge 9 = 9 , with 9q some constant, then this gauge-fixing condition 
remains satisfied throughout the random walk. 

Now consider the gauge theory. Eq. ( |4.4| ) describes a motion in the space of field con- 
figurations that is analogous to that described by the projected equation ( |4.17|) in the toy 



model. In terms of the natural metric ( |A1| ), on the space of gauge configurations, one can 
easily check that the motion is always along directions perpendicular to deformations gen- 
erated by gauge transformations (This is equivalent to satisfying the condition D • A = 0.) 
The question we want to ask is whether the second property of our toy model still holds, 
i.e., is the motion confined to a slice in configuration space? It might be surprising that the 
answer to this question is negative, and one can travel throughout the whole configuration 
space even when restricted to trajectories whose tangents, at every point, are perpendicular 
to the intersecting gauge orbit. This negative answer is perhaps less surprising if one notices 
that there is no obvious gauge-fixing condition similar to 9 = 9$ that is conserved during the 
transverse-projected random walk (|4.4[). Thus, in the gauge theory, there is no equivalence 



of the parameter 9. 

This can be seen most directly by the explicit construction of a trajectory that remains 
perpendicular to gauge transformations at all times, but nevertheless connects two distinct 
points on the same gauge orbit. The trajectory we are going to construct starts at A = 
and remains small all the time, so we can use perturbation theory. Let us denote the small 
parameter by e. Consider first the following trajectory, 



iC,-(x), 0<t<e 



4 u ^ _ ) eC H x ) + (*-e) A(x) , e < t < 2e; , > 

' A ' ' (3e-t)Ci(x) + eD<(x) , 2e < t < 3e; 1 > 

(4e-t)A(x), 3e<t<4e. 

Provided that Cj and Di are transverse, <9jCj = diDi = 0, it is trivial to check that D • A = 
to leading order in e. This means the trajectory is everywhere perpendicular (to within 
0(e 2 )) to the gauge orbits it passes through. 

This trajectory may seem uninteresting, since it is a closed loop that starts at A = 
and ends at the same point. The interest arises when we go to next-to- leading order in e. 
At next-to-leading order, the trajectory ( |B1| ) does not satisfy the condition D • A = 0. For 
example, when e < t < 2e, D • A = e[Cj,Dj]. To correct for this deviation, we need to 
modify Ai in the following way: 

Ai(t, x) = eCi(x) + (t - e)A(x) + (t - e)ai(x), e < t < 2e , (B2) 



where = 0(e), so the term involving aij is of higher order than the other terms. Eq. ( p2| 
satisfies the condition D ■ A = if one places the following constraint on af. 
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didii + e[Ci, Di] = 



(B3) 



This condition can be satisfied by choosing 

a i = -ed i V- 2 [C j ,D j ]. (B4) 

In this manner, one can modify the whole trajectory ( |B1|) so that the condition D • A = 
is satisfied through next-to-leading order. The result is 

(tCt, 

A u^-i ^ + A - e(f-e)$V- a [C,, Dj] , 
^^xj (3 € -t)C i + eD i -(e l + e(t-2e))d i V- 2 [C j ,D j ], 
[ (Ae-t)D l -2e 2 d l V~ 2 [C j ,D 1 ], 

One sees that the trajectory now starts at A; = at t = and runs to Ai = 
—2e 2 d i 'V^ 2 [Cj,Dj} at t = 4e. [Including still higher-order corrections would only change 
this by 0(e 3 ).] At the end point, A is a pure gauge, and, for a general choice of Cj and 
Di, nonzero. Therefore, this trajectory presents a simple example of how, following a trans- 
verse projected trajectory, one can go from one field configuration to a field configuration 
that is gauge equivalent to it. From this result one may show (provided the gauge group 
is semi-simple) that any field configuration in the vicinity of A = is accessible to the 
transverse-projected random walk (|4.4| ). Hence, there can be no analog of the toy model 
"slice parameter" 9 for transverse-projected dynamics in non-Abelian gauge theories. 



< t < e; 
e < t < 2e; 
2e < t < 3e; 
3e < t < 4e. 



(B5) 
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